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CALCULATIONS OF TOTAL CROSS SECTIONS 


FOR CHARGE EXCHANGE IN MOLECULAR COLLISIONS 
INTRODUCTION 

Research accomplished during the first six months of the grant may 
be divided into two main areas. The first concerned the theoretical in- 
vestigation of low energy collisions of nitrogen ions with nitrogen mole- 
cules, Related experimental research was performed at Langley Research 
Center by Dr. A1 Smith of the Instruments Research Division. The effort 
under this grant was to apply theories reported in the literature to this 
specific collision. The calculations aiuL.results are given in the attached 
reports oh pages 2 through 20. • ■ 

The second project was a preliminary search of the literature • on 
collisions of low energy ions wit*h surfaces for the purpose of extending 
the present research. The initial' listing of some reference materials . 
(books, articles, papers, etc.) is attached on page 21. In. addition, a t 
partial literature search was made through the recently acquired Xavier 
Information Program, a computer search of .selected data banks. It is planned 
to use the NASA library search service and to repeat the Xavier search 
when access to more abstracts and data banks, becomes available. 
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MULTISTATE IMPACT PARAMETER METHOD 


This outline of the multistate .impact parameter method 

for the theoretical treatment of molecular charge transfer 

follows closely that given by Flannery, Cosby, and Moran^ 

2 

and uses their notation. The method attempts to take into 
account the vibrational states of the incident ions and of 
the products, x-^ith the target neutral molecules assumed to be 
initially in the ground state. Total cross sections for 
on N 2 are obtained. 

The charge-transfer process -in diatomic nitrogen is 
resonant and symmetric so that the wavefunction "^(t) may 
be expanded in terms of the gerade and ungerade eigenfunctions 


for each electronic state of the quasimoleeular 
complex formed. The atomic separation in the molecule is "f** 

_s» 

and the distance between the tx>70 centers of mass is R. The 
relative phase changes with R because of differences in 


the eigenenergies 


^ (R) . Charge transfer does not occur 

6 , U 


because of an electronic transition but because of this phase- 
3 - 


cnange. 


Assume a straight line trajectory 

A ir 


R = ^ 


a) 





■^.5. 
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where ^ is the impact parameter and v is the incident speed 
along the z-axis . The probability for symmetrical resonant 
charge transfer is 

(f) = h(r'^ ( 2 ) 

for a two-state gerade-ungerade treatment, where 


The charge- transfer cross section is 

r 

cr = H'ty J f (4) 

In order to include nonresonant channels, the wavef unction 

2 r,«i- >s. 

is expanded in terms of ^p^(r) , the molecular eigenfunctions 
of the unperturbed Hamiltonian for- the isolated molecular 
systems. The ^nCt) have eleatrenic, vibrational, and ro- 
tational parts, an^d the eig.enenergies at infinite center-of- 

oL 

mass separation are Eq. 




= Z 

o/ = B,X 


z 

K. 


OL Ci) ^ 




(5) 


The index bd refers to either the direct channels D (no charge 

transfer') or to the exchange channels X. This expansion is 

substituted into the time-dependent Schrodinger equation for 

the internal coordinates. The charge-transfer cross section 
2 

then derived is 

<r = a-ir . j" 1 cj 

o 
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X 

to lowest order. The transition amplitudes C^C^jt) are 
the solutions to the set of ' first-order coupled differential 


eauations 


^ . VCR^ X Cf.« e ^2-. 1,2,.. N 

^T' ifVlij 


where = X when «=^ = D and vice versa. Only the direct 

channel i is occupied initially; 

(8) 

^ = o 

The matrix elements are 

where *z^^) is the vibrational overlap of the incident 

ion in the -z)i vibratxijn^l sta.te and the target neutral in - 

the i) vibrational state. The .rotational motion is neglected 
rx 

since the energies are such that the collision is finished 
before the molecular complex can rotate significantly.^ The 
interaction potential V(R) is assumed to be spherically 
symmetric and can be obtained from the gerade-ungerade splitting 
of the vibrational motion. 

\iR') ^ I f Cr") - (10) 

Then the exchange cross section becomes 
O- = Z-TT ^ ^ 
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The gerade potential is given by the Morse function‘d 




_ -2.^3 CR- Re') -/3^R-Re)- 

^ - 2e, ' ] 


( 12 ) 


where Rg is the internuclear distance corresponding to the 
potential minimum of depth D^, and ^ describes the shape 
of the potential. The modified Lennard- Jones potential^’ ^ 

^ ^ L ( - (. R ) ] 2 R'’ (13) 

where <=< is the polarizability and € and d are Lennard- 
Jones parameters which are kno;«i for the nitrogen potential,^ ' 
gives the interaction energy for the gerade state for this 
reaction. Fitting the Morse curvai_to the Lennard- Jones poten- 
tial allows determination of the parameters Rg = 3.5262 
Dq - 0.06917 eV, and 5903 for nitrogen.^ 

O 

Sato proposed that the ungerade potential may be written 

Si 3 

-V- 2 e. J (2.4) 


The interaction potential V(R) from Eq. (10) is therefore 
o brained to be 


, _ -Z(3(R-0 -fiCR-Re") 

VCR):. <' ^ (15) 


0 
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V(R) is set equal to zero for small -R < R'', where 

R"' = Rg - In 6/|(3 - 2,3995 ^ for N^, since negative values 

of the interaction potential are unphysical. ‘Since D_, 

si- 

and p have been determined, the interaction potential is 
now known. 

Vibrational overlaps are needed in order to obtain the 

transition amplitudes for the charge- transfer cross 

section in Eq. (8) . Overlaps for diatomic ions and neutral- 

molecules in the ground state are given by Nicholls.^ 

Flannery et al^ computed overlaps for excited states of the 

neutral using the Franck-Condon factors^ with -a Morse anhar- 

9 10- 

monic oscillator wa ve function. ’ The excitation defects of. 
various reaction channels were computed from spectroscopic 

4 11 1 9 

energy levels ’ ’ of the' ionic ;and neutral species. They 

found^ that the excitation defects occurred in groups or 
’’bands** separated by approximately one vibrational quantum 
of energy, about 0.25 eV.^^ Therefore the channels could be 
considered to be degenerate within a given band, greatly 
s.im.pl frying the calculations involved. The transition prob- 
abilitries - |c^| were obtained as a function of impact 

parameter ^ by solving the multistate Equations (7) numeri- 
cally by the Adams -Moulton method. The charge transfer cross 
sections cr were obtained from Eq. (6) using these transition 
probabilities by numerical integration over impact parameter 
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using Siiapson’-s Rule with a built-in accuracy parameter. 

The low-velocity approximation can be used for colli- 
sions in which the incident velocity is small enough so that 
the excitation defect is approximately zero. Then 

Equations (7) for the transition probability become 

^ = N/CR'> Z. ( 16 ) 

where B is the number of degenerate final states which have 
zero energy defect with the incident energies. The exchange 
cross section then becomes 

® ^ r ~ -1 

o- = 2-ir X X. \ P J ( 17 ) 

- ■ . <‘fv 

where the are the eigenvalues- and the the asso- 

ciated eigenfunctions of the matrix' which diagonalizes the 

matrix P formed from the P. ; 

sL m 

TT = X-'PX=XPX 08) 

Some representative results of both the low-velocity 
approximat.ion and the full multistate treatment as well as 
some additional experimental data from the literature are 
presented in. the accompanying figures, for comparison with 
•the Xangley data . 
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FIGURE DISCUSSIONS 

Figs. 1-3 present total charge- transfer cross sections 
for i ?2 on N£ obtained at different pressures in the charge 
transfer cell as a function of the square root of the incident 
ion kinetic energy. Each figure contains data points obtained 
at pressures in the CTC of 0.99, 1.33, 1.66, and 1.987 x 10”'^ 
Xorr. Fig, 1 shows cross "sections obtained for ionizing elec- 
tron energy of 16 eV, Fig. 2 for electron energy 17.5 eV, and 
Fig. 3 for electron energy 19.5 eV. 

Figs. 4-7 shox-j the same data of total charge -transfer cross 
sections for H 2 on N 2 as a function of the square; root of 
incident ion kinetic energy. .Each plot contains data for the 
three different ionizing^.electron. .energies (16.0, 17.5, and 
19.5 eV) x^ith a separate figure for each different pressure in 
the CTC. 

Fig. 8 compares the Langley experimental results x^ith 
similar data from the literature and with the theoretical cal- 
culations of Flannery et al.^ Since most of the data for com- 
parison were for ionizing electron energies much higher than 
the present values, the data obtained using the largest elec- 
tron energy of 19.5 eV at an intermediate CTC pressure of 
X.66 xs 10"*^ Torr.x^as 
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chosen for comparison. It can be seen from the previous seven 

figures that the other cross section data falls in the same 

region. The Langley data is shown as filled circles in Fig. 

8. The full multistate treatment (M) and the low-velocity 

approximation (LV) were calculated by Flannery et al^ using 

the distribution of vibrational levels in the -incident ion beam 

corresponding to an electron excitation energy of 22 eV. The 

experimental data obtained by this group is not included in 

Eig. 8 since results were published only for ion energies of 

733 eV and higher. Ghosh and Sheridan^ (GS) state that the 

N 2 ions are in the ground state when leaving their Heil 

type ion source. Pressure in their collision chamber was 

measured by a Knudsen gauge to be._L. 7 x 10 Torr. Potter-^ 

(?) also claims to have ‘Cfely ground state ^2 (X ions in 

2 

his incident beam, although Ghosh and Sheridan disagree with 
his statement. The electron energy in the ion source of 
Stebbln.gs, Turner, and Smith^ (STS) was 200 eV. Utterback 
and Miller^ (UM) used only 22 eV electrons in their ion source, 
and the pressure in their collision region was 2 x 10*“^ Torr. 
Some data ox Berry (B) and Gustafson and Lindholm' (GL) are 

O 

also presented in Fig. 8 for comparison. Nichols and Nitteborn 
(N^a’) measured total cross sections at very low ion energies 
of less than 17 eV. Electron energies in their ion source 
were 95 eV, and their reaction chamber pressure was between 
8 X. 10“^ and 10~^ Torr. 
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Moran et al^ used the multistate impact parameter method 
to obtain cross sections for individual channels and integral 
cross sections summed over all product channels. Theoretical 
cross sections were obtained as a function of the square root 
of incident ion kinetic energy for nitrogen when both the 
incident ion and target neutral were in the ground electronic 
and vibrational state and the resultant product neutral 
and ion in vibrational states equal to (0,1), 

(1,0), and (0,2), and for total vibrational excitation of the 
products. These cross sections presented in Fig. 9 of Ref. 

9 were less than 4 for all cases.,, decreasing to less than 
1 f?‘ for incident ion energies of about 900 eV. They also 
obtained theoretical cross sections-, for total vibrational ex- 
citation with the incident 'Ion in the ground electronic state 
N 2 *^(X~Xg) • and vibrational states 0,1, 2,3 with the target 
neutral in the ground electronic and vibrational state 
N 2 (X^,"=: 0 ) as a function of incident ion kinetic energy. 

•These cross sections presented in Fig. 10 of Ref. 9 were less 
than 10 for incident energies larger than about 25 eV. Since 
-the cross sections x-7ere so much smaller than the other data, 
they were not included in Fig. 8. 
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COLLISION TIMES 


Electronic, vibrational, and rotational transitions will 
be important only when such transitions can take place during 
the time of the collision. The vibrational and rotational 
periods of N 2 are ^ 2 x 10“^^ sec and 2 x 10"^^ sec 

The collision time, i.e.,'the lifetime of the ion-molecule 
complex, is ^ 2 ,4(R/v) 10“-^ sec in atomic units. ^ 

The range of incident energies, 9 eV E 441 eV, or speeds 
5 X 10“^ -d V -d 4 X 10"^ (in atomic units) , is the same order 
of magnitude as the speed range of 7 x l0“'^-«d v-d 6 x 10“” 
given by Flannery et al.^ These speeds involve large inter- 
action radii 10 > R > 5 , for hhe collision. Then the collision 
times are in the range 4.7 x 10 sec *> T' ^ 3.4 x 10*”^^ sec 
Therefore the collision times are much much less than the time 
of rotation, and the collision is finished before the ion-molecule 
complex can rotate. However, the collision time is the same 
order of magnitude as the vibrational time, and a vibrational 
transition is highly likely during the collision. 


1. M.R. Flannery, P.C. Cosby, and T.F. Moran, J, Chem. Phys. 
55^5494 (1973) 



DECAY OF EXCITED STATES OF 

For the measurement of total charge- transfer cross 

sections, consideration should be given to the lifetimes for 

the existence of the various excited states of N 2 . The length 

of the Langley charge transfer cell (CTC) is 1.985 inches, or 

5.0419 cm. For an ion of incident energy 9 eV, or incident 

5 

speed 7.87 x 10 cm/sec, the time spent in the CTC is on the 
order of 6.41 ^^^sec; for an ion of' incident energy 441 eV, 
incident speed 5.51 x 10^ cm/sec, the time is of order 
0.915 JXs&c. . Since the electron ..energies to produce the N 2 
were selected to be bet^- 7 een 16 and 20 eV in the Langley exper- 
iment, only the ground state )Sr 9 ‘ and the first two excited 

states, N 2 '^(A^ 77 ^) and need to be considered. The 

lifetimes of the B— >X radiative transitions are on the order 
of 60 nsec,^’^ so that all the N 9 ' ions formed in the 
electron impact ionization undergo transitions to the ground 
state N 9 (X""^^) within the source region. The lifetimes of the 

"" O 

A — ^ X radiative transitions are on the order of 10 yitsec.^’^’^ 
These ions undergo transitions which occur exponentially with 
distance as where “V is given in the following table 

for the first few vibrational levels.^ The presence of some 
undecayed Ng' (A TTy) states v 7 as included by Flannery et al° 
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since the ions were contained for 9.6 - 14.2 yWsec in their 
time-of-flight apparatus. 

The percent of (electronic and vibrational) ground state 
No (X £o^,'3’=0) in the reactant ion beam as a function of ion- 

£m O 

izing electron energy from 0 to 100 eV is presented in Fig. 4 
of Ref. 6. For energies of 16, 17.5, and 19.5 eV (those used 
at Langley) , it can be deduced from this figure that between 
80 and 90 percent of the incident ion beam must be in the 
ground state 'N2*^(X^2j, ^*’=0) . 

O 

i) T 

0 15.5 

1 13.7 

2 12.2 

3 ’ 11,0 ■ 

4 9.9 

5 9.1 

6 8.4 
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